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Abstract 

The operation of replacing every vertex of an r-regular lattice H by a complete graph 
of order r is called clique-inserting, and the resulting lattice is called the cliquc-insertcd- 
lattice of H. For any given r-regular lattice, applying this operation iteratively, an infinite 
family of r-regular lattices is generated. Some interesting lattices including the 3-12-12 lattice 
can be constructed this way. In this paper, we reveal the relationship between the energy 
and resistance distance of an r-regular lattice and that of its clique-inserted-lattice. As an 
application, the asymptotic energy per vertex and average resistance distance of the 3-12-12 
and 3-6-24 lattices are computed. We also give formulae expressing the numbers of spanning 
trees and dimers of the A:-th iterated clique-inserted lattices in terms of that of the original 
lattice. Moreover, we show that new families of expander graphs can be constructed from 
the known ones by clique-inserting. 

Keywords: Clique-inserting, Lattice, Kirchhoff index, Dimer, Spanning Tree, Expander 
graph 



1 Introduction 

The study of 2-dimensional lattice models in statistical physics has a long history. One kind of 
lattices that have been paid a lot of attention to are constructed by replacing each vertex of r- 
regular lattices with a complete graph of order r.(To avoid triviality, we assume r > 3 throughout 



the paper.) Such lattices include the martini [18j[2TJ[23] , the 3-12-12 [19j[2TJ[22] , the 3-6-24 [Tl] 
and the modified bath room lattices [21] . Following [29] , the operation of transforming each 
vertex of an r-regular graph to an r-clique (complete graph of order r) is called clique-inserting, 
and the graph obtained this way is called the clique-inserted-graph of the original graph. From a 
given r-lattice, the operation of clique-inserting can also be performed, and the resulting lattice 
is called the clique-inserted-lattice of the original lattice. 

Throughout this paper, we always assume that G denotes an undirected simple graph. Note 
that in the language of graph theory, the clique-inserting operation on a graph G can be described 
as taking the line graph of the subdivision graph of G. For any given regular lattice, by iterating 
this operation, a set of hierarchical regular lattices, namely, iterated clique-inserted-lattices can 
be obtained. For instance, based on the hexagonal lattice, the 3-12-12, 3-6-24 and 3-6-12-48 
lattices (refer to |11| for definitions of these lattices) can be generated by clique-inserting. In 
this case, the clique-inserting operation on a lattice is equivalent to a star-triangle transformation 



on the subdivision lattice of the original lattice. By this observation we compute some physical 
and chemical indices of the 3-12-12 and 3-6-24 lattices from the hexagonal lattice which has 
been well studied. This approach can also be applied on regular lattices of degree greater than 
or equal to four, such as the modified bath room lattice (the clique-inserted-graph of square 
lattice). 

For a regular lattice H, let C(H) denote the clique-inserted-lattice of H. Denote the sequence 
of lattices obtained by iterating clique-inserting operation as follows, C°(H) = H, C l (H) = 
C(H), C 2 {H) = C(C(H)), C k {H)=C{C k - l {H)), ... where k = 1,2,.... In this paper, we 
consider the lattices produced by the operation of clique-inserting on regular lattices with free, 
cylindrical and toroidal boundaries. We will discuss the energy per vertex, average resistance 
(Kirchhoff index over the number of edges) and the entropy of dimer and spanning trees of such 
lattices. We will also use the operation of clique-inserting to construct new families of expander 
graphs from known ones. 

The dimer problem of lattice has been attracting the attention of many physicists as well as 
mathematicians. For some classical works, we refer to (6}[20}[22). Cayley [I] and Kirchhoff [13] 



presented the problem of enumeration of spanning trees of graphs, and further work in statistical 
physics has been appearing in both physics and mathematics literature. For a good survey, the 



reader is referred to 19 . The study of resistance distance is initiated by Klein and Randic 14 



and the related index named "Kirchhoff index" was well studied in 110,24 . Gutman |8| defined 



the energy of a graph G as the sum of the absolute values of eigenvalues of G |9j. Yan and 



Zhang 27 proposed the energy per vertex problem for lattice systems and showed that the 
energy per vertex of 2-dimensional lattices is independent of various boundary conditions. For a 
comprehensive survey of results and common proof methods obtained on graph energy, see [15] 
and references cited therein. Expander graphs were first defined by Bassalygo and Pinsker in 
the early 70's. They are regular sparse graphs with strong connectivity properties, measured 



by vertex, edge or spectral expansion as described in 12 . The property of being such graphs 
turns out to be significant in various disciplines. For example, expansion is closely related to 
the convergence rates of Markov Chains and so it plays a crucial role in the study of Monte- 



Carlo algorithms in statistical mechanics 12 . Due to the usefulness of expander graphs in 
complexity theory, computer networks and statistical mechanics, their constructions have been 
studied extensively. 

The rest of the paper is organized as follows. In section 2 and 3, the expression of the energy 
and Kirchhoff index of A;-th iterated clique-inserted-lattices of regular lattices are discussed, 
respectively. As an application, we compute the energy per vertex and average Kirchhoff index 
of the 3-12-12 and 3-6-24 lattices. In section 4, we show that, given Zh as the entropy of spanning 
trees of a r-regular lattice H, the entropy of spanning trees of C k (H) (the k-th iterated clique- 
inserted graph of H) is given by r~ k {zu + s In r(r + 2)) where s = (r/2 — l)(r k — l)/(r — 1). We 
will also show that when H is cubic, the free energy per dimer of the C k (H) is ^ In 2. In section 
5, we will construct new families of expander graphs from known ones by clique-inserting. 

2 Asymptotic Energy 

Let G = (V(G), E(G)) be a graph with vertex set V(G) = {t>i,t>2, • • • ,v n } and edge set E(G). 
The adjacent matrix of G, denoted by A(G), is an n x n symmetric matrix such that ay = 1 



2 



if vertices Uj and Vj are adjacent and otherwise. Line graph L(G) of G, is a graph such that 
each vertex of L(G) represents an edge of G and two vertices of L{G) are adjacent if and only if 
their corresponding edges of G share a common end vertex in G. The subdivision graph S(G) 
of a graph G is the graph obtained by inserting a new vertex into every edge of G. It is easy to 
see that C(G) = L(S(G)). The energy of a graph G with n vertices, denoted by £(G), is defined 
by £ (G) = J27=i l^i(^)l) where Xi(G) are the eigenvalues of the adjacency matrix of G. The 



asymptotic energy per vertex of G 



27 



is defined by lim 



g(G) 



\V{G)\Uoo WW 



Lemma 2.1. '2V Suppose that {G n } is a sequence of finite simple graphs with bounded average 
degree such that lim |V(G n )| = oo and lim /(n n M = h / 0. If {G'\ is a sequence of spanning 

subgraphs of {G n } such that lim n \-Jfp^ ^ — ^ = 1, then lim r§7^ri = h. That is, 

n— >oo \ v \S"n>)\ n— >oo \ v \'- r n)\ 

G n and G' n have the same asymptotic energy. 



Lemma 2.2. 291 Let G be an r -regular graph with n vertices and m edges. Suppose that the 
eigenvalues of G are X\ = r > A2 > ... > A n . Then the eigenvalues of the clique-inserted graph 

C{G) of G are r 2 ±-\A^+ 4 ( A ' +1 ) ^ ^_ 1^2,..., n, besides —2 and each with multiplicity m — n. 

From Lemma 2.2, we immediately obtain the following corollary. 

Corollary 2.3. Let G be an r-regular (r > 3) graph with n vertices and eigenvalues Ai > A2 > 
... > A n , the energy of the clique-inserted graph of G is 

£{C(G)) _ £ r-2 + ^ + 4(A, + l) + £ y^ + 4(A, + l)- r + 2 + „ 

1=1 8=1 

We will use this result to calculate the asymptotic energy per vertex of the 3-12-12 lattice 
and its clique-inserted-lattice in the rest of this section. 



2.1 3-12-12 lattice 



Our notation of hexagonal lattices follows 27 . The hexagonal lattices with toroidal condi- 
tion, denoted by H (n, m), are illustrated in Figure 1, where (a\, a*), (02, a^), ■ ■ ■ , (a m +i, a* m+1 ), 
{bi,b\),{b 2 ,bl),...,(b n+1 ,b* n+l ) are edges in ff*(n,m). 




Fig. 1: H*(n, m) with toroidal boundary condition 
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By the definition of clique- inserted-lattice, it is easy to see that each 3-12-12 lattice with 
toroidal boundary condition is a clique-inserted- graph of H (n, m), denote as T l (n,m) (see 
Figure 2(a)). Note that (<n, aj), (a 2 , aj), . . . , (a m+ i, a^ +1 ), (&i, 6J), (6 2 , b* 2/ 



, . . . , (6„+i, 



arc 



edges in T*(m, n). If we delete edges (61, 6*), (62, ^1), • • • 5 (b n +ii ^n+i) from T*(n, m), then the 
3-12-12 lattice with cylindrical boundary condition, denoted by T c (n,m) (see Figure 2(b)) can 
be obtained. If we delete the edges (ai,o*), (02, a|), . . . , {a m+ i,a* m+l ) from T c (m, n), then the 
3-12-12 lattice with free boundary condition, denoted by T^(m,n) (see Figure 2(c)) can be 
obtained. 




Fig. 2: The 3-12-12 lattice T*(n,m) (left), T c (n,m) (middle), and T f (n,m). 



Theorem 2.4. For i/ie 3-12-12 lattices T t (n,m),T c (n,m), and T^{n,m) with toroidal, cylin- 
drical and free boundary conditions, 



lim 

n,m— » 00 



£(T*(n,m)) 
6mn 



lim 

n,m— > 00 



£{T c {n,m))) 



6mn 

r-2n r 2n 



lim 

n,m— » 00 



£(Tf(n,m)) 



3 + 24^2 / y ( V 13 + 4 V 73 + 2 cos x + 2 cos y + 2 cos(x + y) 



'0 JO 

+^/l3 - 4 ^3 + 2 cos x + 2 cos y + 2 cos(x + y))d x d y « 1.4825 

i/ia£ is, the 3-12-12 lattices T t {n,m),T c {n,m), andTf(n,m) with toroidal, cylindrical, and free 
boundary conditions have the same asymptotic energy (~ 8.895mn). 

Proof. Note that almost all vertices of T c (m, n) and T*{m, n) are of degree 3. Since T* (m, n) 
and T c (m, n) are spanning subgraphs of T*(m, n), by Lemma 2.1 we have 



lim <S ( T ( n ' m ))) £(T c (n,m))) 

n,m— >oo Qnin n,m—>oo 



6mn 



lim 

n,m— > 00 



£(T'(ra,m))) 



6mn 



It is shown in 27 that the eigenvalues of H (n,m) are: 



± y 3 + 2oo8^ + 2coB^ [ + 2coB^ + ^ [ j,0<i<n 1 0<i< 
Since T*(n, m) is the clique-inserted graph of H l {n, m) , we have 



2&JL 



2in , 2j7r 



m. 



£(r'(n,m))) 



EE\ 

7=0 \ 



i=0 



/ 277T 277T / 2Z7T 277T 

13 + 4W3 + 2cos + 2 cos J — +2 cos ( — — + ' 



n + 1 



m + 1 



n + 1 m + 1 



EE 

i=o j=o 



1 



/ 2i7T 277T / 2Z7T 2771 . 

13 - 4W3 + 2 cos + 2 cos J — + 2 cos H —— + 2mn. 



n + l 



m + 1 



77+ 1 777 + 1 
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Thus, the average energy per vertex of 3-12-12 lattice can be expressed as 

£(T l (n m)) 1 1 f 2w f 2w I ~i^^^^^^^^ = ^^^^^^^^ = - 

lim v / ' n = - + — / / (a/13 + 4 a/3 + 2cosx + 2cosy + 2cos(x + 
n,m->oo 6mn 3 24:7t z J J v 



V) 



+\J 13 - 4 a/3 + 2 cos x + 2 cos y + 2 cos(x + y))d x d y « 1.4825 
and we complete the proof of the theorem. □ 



2.2 3-6-24 lattice 




Fig. 3: The 3-6-24 lattice S t (n,m) (left), S c (n,m) (middle), and (n,m) (right). 



The clique-inserted-graph of T t (m,n) is a lattice with toroidal boundary condition, denoted 
by 5*(m, n), illustrated in Figure 3. Note that (ai, a*), (02, aJj), . . . , (a m+ i, a* t+1 ), (61, (62, W>); 
. . . , (b n+ i, 6* +1 ) are edges in 5*(m, n). If we delete edges (61, 6*), (b 2 , b* 2 ), (6 n+ i, from 
S**(n, m), then the 3-6-24 lattice with cylindrical boundary condition, denoted by S c (n,m) 
(see Figure 3(b)) can be obtained. If we delete edges (a±, a*), (02, a 2 )j ■ ■ ■ , (a m +i, a* n+ i) from 
S c (m,n), then the 3-12-12 lattice with free boundary condition, denoted by S^{m,n) (see Fig- 
ure 3(c)) can be obtained. 

Theorem 2.5. For 3-6-24 lattices S^m, n), S c (m,n), and S^(m,n) with toroidal, cylindrical 
and free boundary conditions, 



lim 



£(S*(m,n)) 



n,m-s>oo 18mn 



lim 



£(S c (n, m)) 



n,m-s>oo 18mn 



lim 

n,m—>- oo 



£(Sf(n,m)) 
18mn 



2 >>27r (>27r / / — 

= T2^J J ( V 15 + 2 V 13 + 4 ^ 3 + 2 cos x + 2 cos y + 2 cos(x + 



V) 



+\/ 15 + 2a/ 13 - 4 a/3 + 2 cos x + 2 cos y + 2 cos(x + y) 



+y 15 - 2a/ 13 + 4 a/3 + 2 cos x + 2 cos y + 2 cos(x + y) 

+-\y 15 — 13 — 4a/3 + 2 cos 2 + 2 cos y + 2 cos(x + y))d x d y + — + + 6 



1.4908 



that is, the lattices S t (n,m)), S c (n,m))), and S^{n,m)) with toroidal, cylindrical, and free 
boundary conditions have the same asymptotic energy (m 26.8344mn/ 
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Proof. Note that almost all vertices of S c (m, n)) and (m, n)) are of degree 3. Since S$ (to, n) 
and S c (m,n) are spanning subgraphs of S t (m,n), by Lemma 2.1 we have 

£(S*(n,m)) ,. £(S c (n,m)) ,. S(SUn,m)) 

hm = hm = hm 

n,m->oo 18mn ri,m->oo 18rnn n,m-5>oo 18rnn 

The energy of the the clique-inserted-graph of 3-12-12 lattice can be obtained by 



f(S*(n,m)) = EEa 
i=0 i=0 \ 



15 + 2 



13 + 44/3 + 4 cos h2cos — h2cos( — — + 



n + 1 



m + 1 



n + 1 m + 1 




13-4W3 + 4cos + 2 cos J -— + 2 cos + 



n + 1 



m + 1 



n + 1 m + 1 



+EE 

i=0 j=0 



, / , 2Z7T 2?7T / 2i7T 27 7T 

13 + 4j3 + 4cos — — + 2 cos - + 2 cos - + 



n + 1 



m + 1 



n + 1 m + 1 



+EE 

i=0 j=0 
+V / 5mn + V4~3mn + 6mn 



=0j=0 \ \ 



/ 2z7r 2?7r / 2iir 2?7r 

13-4W3 + 4cos h2cos — — + 2cos + — — 

V n + 1 m + 1 \ n + l m + 1 



The average energy per vertex of the clique-inserted-lattice of the 3-12-12 lattice can be obtained 
by 



lim 



S(S\n,m)) 



J J ( y 15 + 2y 13 + 4^3 + 2 cos x + 2 cos y + 2 cos(x + yj 



n,m-s>oo 18mn 



72vr 2 



+\/ 15 + 2a/13 - 4^3 + 2 cos x + 2 cos y + 2 cos(x + y) 



+Y 15 - 2\j 13 + 4^3 + 2 cos x + 2 cos y + 2 cos(x + y) 

+y 15 - 2 Y^13 - 4^3 + 2 cos x + 2 cos y + 2 cos(x + y))d :E dj / + — + ^ + 6 



1.4908 



and we finish the proof of the theorem. 



□ 



3 Average Resistance 

Kirchhoff index Kf(G) of a graph G is defined as the sum of resistance distance between all 
pairs of vertices of G. That is, Kf(G) = R(u, v), where R(u, v) denotes the resistance 

{u,v}GV(G) 

distance between vertices u and v of graph G. Let Kf{C(G)) = f^-Kf(G) denote the average 
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Kirchhoff index, that is, the average resistance distance between all pairs of vertices of G. Unlike 
the case of energy per vertex, an edge deletion of the graph may change the average resistance 
dramatically. For example, if we delete a cut edge of the graph, the resistance distance will 
become infinity. Thus we only consider lattices with toriodal boundary condition here. 

Lemma 3.1. ^Let G be a connected r -regular graph with n > 2 vertices. Then 

Kf(L(G)) = r -Kf{G)+ {r ~^ )n2 . 
Lemma 3.2. Let G be a connected r-regular graph with n > 2 vertices. Then 

A7(s(G))=^A7( G) + (r2 - 4) ; 2 + 4 " . 

Combine above two lemmas, the following results are straightforward. 
Proposition 3.3. Let G be a connected r-regular graph with n > 2 vertices. Then 

Kf(C(G)) = r(r + 2)2 A 7 (G) + r3n2 " 2rn2 16 +4rn " 4n2 
Corollary 3.4. Let G be a connected 3-regular graph with n>2 vertices. Then 

Kf(C(G))^Kf(G) + ^±^ 

3.1 3-12-12 lattice 

Theorem 3.5. For the 3-12-12 lattice T*(m, n) with toroidal boundary condition, 
Kf(T\m,n)) « 204.9788 + ^ + ^ + % + ^ + + ^ , 



lim Kf{T\m,n)) = 0.2361111.... 

n,m— >oo 

Proof. It is shown in [28] that for the hexagonal lattice H (m,n) with toroidal boundary 
condition, 

Kf{H t {m, n)) « 10.9322(m + l) 2 (n + l) 2 . 
Therefore, by Corollary 3.4, we have 

„ t , » 17(m+ l) 2 (ra+ l) 2 +6(m+ l)(n+ 1) 

Kf(T l {m, n)) = 204.9788 + — ^ — ; K ; v — A ; 

Thus, the average Kirchhoff index of T*(m, n) is given as follows, 

lim Kf(T\m,n)) = lim ^g^^ = 0.2361111... 

□ 
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3.2 3-6-24 lattice 

Theorem 3.6. For the 3-6-24 lattice S t (m,n) with toroidal boundary condition, 

KHS'^.n)) p, 3843.3525 + 1887 ( m + 1 ) 2 (" + !) 2 + 522(m + l)(n + 1) ^ 

16 

W](S\m,n)) = 0.7280093.... 
Proof. By Corollary 3.4 and Theorem 3.5, we have 

K /(SV, „)) p, 3843.3525 + lB87( m + 1)> + 522(,n + !)(„ + 1) ^ 

Thus, the average Kirchhoff index of S t {m, n) is given as follows, 

lim Kf(S\m,n)) = lim ^^^11 = 0.7280093.... 

ra,m->oo J v v ' /y n.m-s-oo /18(m+l)(n+l)^ 

□ 

4 Spanning Tree and Dimer 
4.1 Spanning Tree 

Let Nst(G) denote the number of spanning trees on G. It is known that Nst(G) has asymptotic 
exponential growth. Define the quantity zq by 

zq = lim — In Nst(G) 

n— >oo n 

This quantity is known as bulk, or the thermodynamic, limit, and we denote zq by zl for lattices. 
The following lemma indicates the relation between the number of spanning trees of a regular 
lattice and that of its fc-th iterated clique- inserted-lattice. 

Lemma 4.1. Jig 5) / Let G be an r-regular graph with n vertices. Then the number of span- 
ning trees of the iterated clique-inserted- graphs C k (G) of G can be expressed by NsT(C k (G)) = 
r ns-k^ r + 2) ns + k N ST (G) where s = (r/2 - l)(r k - l)/(r - 1). 

Therefore, we obtain the following result. 

Proposition 4.2. Let H be an r-regular lattice. For C k (H)(k = 0, 1,2, ..), the rate of growth 
of number of spanning trees, z C k^ is given by r~ k (zn + slnr(r + 2)) where s = (r/2 — l)(r k — 
l)/(r — 1) and zh denotes the rate of growth of spanning trees of H . 

The next Theorem implies that the boundary condition does not affect the bulk limit of a 
lattice. 
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Theorem 4.3. fFq Let (G n ) be a tight sequence of finite connected graphs with bounded average 
degree such that 

lim | V{G) n T 1 | {x e V(G' n y,deg G , n (x) = deg Gn (x)} \= 1, 

n— too 7 

then lim | V(G) n \~ l log t(G') = h. 

n— >oa 

For hexagonal lattice, Zh c is 0.8076649... as shown in |19| . Thus, by Theorem 4.2 and 
Theorem 4.3, we have that for 3-12-12 lattice and the 3-6-24 lattice with toroidal, cylindrical 
and free boundary condition, 

23-12-12 = 0.7205633... 
■Z3-6- 24 = 0.69 1 52 95.... 

4.2 Dimer Problem 

Let M(G) denote the number of dimers (perfect matchings) of G. The free energy per dimer of 
G, denoted by Z G , is defined as Z G = lim f In M(G). Given the number of vertices and edges 

n— >oo 

of a connected graph, the number of dimers of the graph and that of its line graph have the 
following relation. 

Lemma 4.4. |5]/ Let G be a 2-connected graph of order n and size m, where m is even. Then 
M{L{G)) > 2 m - n+1 , where the equality holds if and only if A(G) < 3. 

With this general result, we can readily obtain the following. 

Theorem 4.5. Let H be a cubic lattice H with toroidal condition, the free energy per dimer of 
C k (H) (k = 1,2,3,..; is I In 2. 

Proof. Let H be a cubic lattice with n vertices, then the edge set of H is of order Note that 

C fe_1 (i^) has 3 k ~ 1 n vertices and edges and the subdivision of C fc_1 has | • 3 fc_1 n vertices 
and 3 k n edges. Since C k (H) is the line graph of the subdivision of C k ~ 1 (H), by Lemma 4.4 we 
have Z rk (G) = lim J- In 2 3 fc n-|-3 fc n+i = i ln2 . □ 

Theorem 4.6. Let i?*(m,n) be the k-th iterated clique-ins erted-lattice of the hexagonal lattice 
LL^im^v^with toroidal boundary. For the the lattices i2*(m,n), R c (m,n) and Rf(m,n) with 
toroidal, cylinder and free boundary conditions, Z^t^^ = Z^c^ m ^ = Z R f^ m ^ = 5 In 2. 

Proof. Note that R c (m,n ) (Rf(m,n)) can be considered as the line graph of a graph which 
differs from S(C k ^ 1 (H t m,n)) by small number(small in the sense that the number is o(mn) as 
m,n approach infinity) of edges. Therefore, the results follows by applying Lemma 4.4. □ 

In general, when a cubic lattice is a line graph, the free energy per dimer of plane lattices are 
the same as that of the corresponding cylindrical and toroidal lattices. However, this may not 
be true when a cubic lattice is not a line graph. The hexagonal lattice is such a counterexample 
as shown in |26|. 
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5 Expander graphs 



Let dc(vi) be the degree of vertex v% of G and D{G) = diag(d,G(vi), da(v2), ■ ■ ■ , da[v n )) be the 
diagonal matrix of vertex degree of G. The Laplacian matrix of graph G is L{G) = D{G) — A{G). 
The eigenvalues of L(G), denoted by Hi < Hi < ■ ■ ■ < [in are called the Laplacian spectrum of 
G. It is well known that \xi is called the spectral gap (or algebraic connectivity) of G, and fi2 is 
greater than if and only if G is a connected graph. Furthermore, (12(G) is bounded above by 
the vertex connectivity of G. Here we use the spectral gap to quantify the expansion property 
that is, a family of regular graphs is a an expander family if and only if it has a positive lower 
bound for the spectral gap, and the larger the bound the better the expansion. This character- 
ization can be formulated to a formal definition as the following. 

An infinite family of regular graphs, G\, G2, G3, . . ., is called a family of e-expander graphs (12] , 
where e > is a fixed constant, if (i) all these graphs are /c-regular for a fixed integer k > 3; (ii) 
(12 > £ for i = 1, 2, 3, . . .; and (iii) rij = |V(Gj)| — > 00 as i — > 00. Note that for regular graphs, 
\ii = r — Xi for i = 1,2, ... ,n. Then Lemma 2.2 implies that 



, (r(r] , r + 2-^(r + 2f-A(i 2 (G)) 
fi2(^{G)) = . 



Denote the function iteration of f{x) = r + 2 V(^+ 2 ) 4rr f 1 ^) = f{x) and f k+1 (x) = f(f k (x)) 
for k = 1, 2, 3, Then we can readily obtain the following result. 

Corollary 5.1. Suppose G\,G2,Gs, . . ., is a family of r -regular e-expander graphs. Then 
C k {G\), C k (G2), C k {Gz), . . . , is a family of r-regular f k (e)- expander graphs. 

Lattices with good expansion property can be used to design the network with optimal syn- 
chronizability and fast random walk spreading [3] . Though the best expansion for combinatorial 



expanders is still unknown 16 , there is an upper bound on the spectral gap. Given by the 
spectral gap of the Bethe lattice, Alon and Boppana derived the asymptotic upper bound as 
follows. 

Theorem 5.2. fToy Let G n be a family of r-regular graphs, where r is fixed and n is the number 
of vertices of the graphs. Then 

lim H2{G n ) <r-2y/r-l. 

n— >oo 

A r-regular graph has fj,2 — r ~ 2\/r — 1 (or equivalently say, A2 < 2y/r — 1) is called a 
Ramanujan graph. By Alon-Boppana Theorem, Ramanujan graphs are a family of expanders 
which achieve the largest possible spectral gap and thus have the best possible expansion from 
a spectral point of view. 

In order to introduce some existing families of Ramanujan graphs and demonstrate our construc- 
tion of new expander families based on that. The following definitions need to be introduced. 

Definition 5.3. Given a group T and a generating set S C T \ {1} such that 5" _1 = S, the 
Cayley graph Cay(T, S) is defined to have vertex set X such that x,y £ V are adjacent if and 
only if x~ 1 y G S. 



10 



Definition 5.4. The group PGL(2,Z q ) is the quotient group of all two by two invertible matrix 
over Z 9 modulo {(q^),A^0} 

Definition 5.5. The group PGL(2,7* q ) is the quotient group of all two by two invertible matrix 
over 1j q with determinant one, modulo {(oi)i('o 1 

Given two distinct odd prime number p and q such that q > 2v2. Define the generating set 
to be 

S = ( a ° + QlX + a3V ~ aiV + a2 + a3X ^\ (5 i) 

\-a\y - a 2 + a 3 x a - a\x - a 3 y J 1 " 

that satisfies the following 

(i) x, y are odd primes satisfying (x 2 + y 2 + 1) = mod q. 

(ii) a 2 , + a\ + a\ + a| = p where a« G Z and ao is a positive odd number if p = 1 mod 4 while if 
p = 3 mod 4 then oo is even and the first non-zero ai is positive. 

Define X p,q to be the Cayley graph Cay(PSL(2,Z q ), S) when p is a square module q and 
Cay(PGL(2,Z q ), S) otherwise. Lubotzky, Phillips and Sarnak provided the above construction 
and proved that X p,q is a family of Ramanujan graphs of degree p + 1 for each prime odd p by 
using some deep result in number theory [4] . Chiu [2] extended this approach to the case of 
p = 2 by considering the generating set of 



(5.2) 



1 0\/2 + V=2 ^26 \/2-7=2 -7=26 
-lM 7=26 2-V=27 1 -V^26 2 + 7=2 



By acting this generating set on the PGL(2,Z q ) or PSL(2,Z q ) as before, a family of cubic 
Ramanujan graphs X 2,q is explicitly constructed. 
It is known that 14], for a fixed real number < 7 < 1/6 and sufficiently large q, the spectral 



Fig. 4: Cubic Ramanujan graph X 2 ' 3 (left) and C(X 2 ' 3 ) (right). 

gap of X p,q is lower bounded by e{r) = {p + 1) — p^ +1 —p&~' 1 . Therefore, for a fixed odd prime 
p, C(X p ' q ) is a p+3 V^(p+ 3 ) 4 ((p+ 1 ) p 5 p b ^ -expander family with degree p + 1. Generally, 
C^X^' 9 ) is a f k ((p + 1) — ps" 1 " 7 — p6 _ i')-expander family. Similarly, as we apply clique-inserting 
iteratively on the expander family X 2 ' q , new families of cubic expander family can be obtained. 
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